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Numerical studies of fractal curves in the plane often focus on subtle geometrical properties such 
as their left passage probability. Schramm-Loewner evolution (SLE) is a mathematical framework 
which makes explicit predictions for such features of curve ensembles. The SLE prediction for the 
left passage probability contains the Gauss hypergeometric function 2 + 1 - To perform computational 
SLE studies it is therefore necessary to have a method for numerical evaluation of 2+1 in the 
relevant parameter regime. In some instances, commercial software provides suitable tools, but 
freely available implementations are rare and are usually unable to handle the parameter ranges 
needed for the left passage probability. We discuss different approaches to overcome this problem 
and also provide a ready-to-use implementation of one conceptually transparent method. 


I. INTRODUCTION 

Quantifying the geometry of fractal boundaries and paths has a long history in the study of phase transitions and 
lattice models. In 2000, Schramm’s work on loop-erased random walks introduced a new formalism devoted to this 
task [1]. Partially due to this formalism, called Schramm-Loewner evolution (SLE), it became possible to prove the 
existence and numerical values, e.g., of the 2D percolation critical exponents [2, 3]. Since SLE describes geometrical 
features of fractal curves, beyond the numerical value of the fractal dimension, it became a useful tool to understand 
several statistical physics models involving random paths and curves [4-7]. Authoritative mathematical reviews of 
SLE can be found, e.g., in Refs. [8-11]. 

Since SLE implies several geometrical properties of the associated curves, it is not clear a priori whether a given 
distribution of fractal curves can in fact be described by SLE. Therefore there have been a number of computer 
simulation studies testing numerically the compatibility of the properties of different ensembles of curves with the 
predictions of SLE [12-18]. The present work is based on a note in Ref. [19]. 

One of the most conclusive and subtle tests of SLE properties deals with the curves’ left passage probability, in a 
certain geometrical setting (known as chordal SLE). For a curve originating in the origin and going to infinity in the 
upper half plane, Schramm proved [20, 21] that the probability P K (</) that the curve passes to the left of the point 
zq = Rq exp(i</>), see Fig. 1, is given by: 

pm) = \ + vct/4 1 ) com 2Fi (+ *+ ~ c ° m2 ) ■ (1) 

Here, k > 0 is the single parameter characterising the geometry of SLE curves [for Eq. (1), 0 < n < 8] and known [9] 
to be related to the curves’ fractal dimension df by 

df = minjl + k/8, 2}. (2) 


As seen from Eq. (2), df = 1 corresponds to k = 0 and df = 2 is recovered for k > 8. We note that the left passage 
probability P K (</>) of an ensemble of curves depends only on n and the angle with the real axis <fi . Figure 2 shows 
P k {4>) as function of </>, for different values of n. For <fi = tt/2, one has P k (tt/2) = 1/2, independent of n, as expected 
by symmetry. For k = 8, the curves are space-filling and Ps(<t>) = 1/2, independent of </> [21]. 

The data in Fig. 2 has been obtained by evaluating Eq. (1) with Mathematica [22]. However, in numerical tests of 
SLE predictions, one needs to evaluate Eq. (1) many times [23-26] and it can be advantageous to do the evaluations 
within the simulation run. In many cases, corresponding simulations are performed using the C++ programming 
language and the associated libraries [27]. The required numerical evaluations of the cotangent and gamma [45] 
functions can be handled by the C++ standard library, the gnu scientific library (GSL) [35] and the Boost libraries 
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FIG. 1. Sketch showing the geometrical definitions involved in the left passage probability formula in Eq. (1). P K (cj>) is the 
probability of a chordal SLE curve with parameter k to pass to the left of a point with argument (/>■ 


[28]. Although GSL and Boost, provide implementations of 2 -F 1 these do not cover the parameter regimes required 
for the left passage probability. 

The remainder of this work is structured as follows. In Sec. II, we describe different approaches of how to numerically 
evaluate 2 -Pi as needed for the left passage probability. Use of argument transformations is reviewed in Sec. IIA and 
integrating 2 -Pi’s integral representation is discussed in Sec. IIB. We draw conclusions in Sec. III. 


II. METHODS 

We note that the evaluation of 2 -P 1 is in general a difficult numerical task requiring a sophisticated combination of 
various techniques [32-34]. Here we focus on methods useful for the parameter regime relevant to the left passage 
probability formula in Eq. (1). 


A. Argument transformations 


Many implementations of 2 -P 1 , such as those in the GSL, Boost, or Cephes [36] libraries only support arguments 
|^| < 1. However, for the left passage probability in Eq. (1), 2 = — cot(</>) 2 , therefore 2 can in general take any non¬ 
positive value. To use such an implementation for [ 2 ] > 1, known relations can be applied to transform the problem 
back to the unit disk [33, 34]. An implementation in a large public library which includes these transformations is 
SciPy [43]. Though the SciPy interface is in python, the underlying algorithm is in c and was adapted from the 
Cephes mathematical library [36] to support z < — 1. 

For — 00 < z < — 1, Ref. [33], proposes to use the following transformation (Abramowitz and Stegun (A&S) [44] eq. 
15.3.8): 


2 F 1 (a, b,c ; z) 


= (!-*) 
+ (!-*) 


—a 


-b 


r(c)r(6 — a) 
r(6)r(c-a) 2 
r(c)T(a — b) 
r(a)T(c — b) 2 


Fi(a,c — b,a — b + 1; (1 


Fi(b, c — a, b — a + 1; (1 




( 3 ) 

( 4 ) 


SciPy, which we consider for the comparison in Fig. 3(b), uses, e.g., for —00 < z < —2, the following relation (A&S 
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FIG. 2. Graph of Eq. (1): For different values of k, Pk(4>) is shown versus the angle <$>/n. The corresponding k values are 
shown next to the curves. Note that, e.g., Pg(</>) = 1/2, Pi{ij>) = 1 — 0/tt, and = 1 + [sin(2</>)/2 — <j>\/n [21]. 


[44] eq. 15.3.7): 


2 F 1 (a,b,c;z) = (-z) — ^\ 2 F 1 (a,l - c +a,l - b + a;l/z) 

1 (b)L {c — a) 

+ 1 - c + 6,1 - O + 6; l/z). 


(5) 

( 6 ) 


When using these transformations care must be taken to avoid special parameter values that are numerically difficult. 
For example, when b — a is an integer (which will happen for k = 8/3) T(& — a) or T(a — b) will diverge. In these cases, 
or when z is close to —1, one could use the transformations A&S 15.3.4 or A&S 15.3.5 [43, 44]. 


B. Integral representation 

As an alternative to using argument transformations, we next apply the method of special function evaluation 
proposed by Schwartz [29-31]. 2 Fi can be expressed as 

2 Fi(a, 5, c; z) = ^ J q t b ~\l - t) c ~ b -\l - zt)~ a dt, (7) 

where Re(c) > Re(6) > 0, |arg(l — z)| < 7r, and T is the gamma function [37-40]. For computing numerically the 
integral in Eq. (7) we use the method proposed in Refs. [30, 31]. The main point of this method is to look at integral 
representations of special functions. Where necessary, they are reformulated as infinite integrals over quickly decaying 
integrands. These integrals are then computed with the trapezoidal rule: 

/ OO __ 

f(x)dx « 'Y^hf{kh), ( 8 ) 

-°° k 

where h denotes the step size. In computational practice, after reaching a desired accuracy, the sum in Eq. (8) is 
truncated (see Ref. [30], and for a recent mathematical treatment see Ref. [41]). This allows to move the singularities, 
for t = 0 and t = 1, of the integrand in Eq. (7) to infinity. Therefore they do not affect the truncated sum of Eq. (8) 
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[31]. For example, this can be achieved by a changing the integration variable (from t to u ) as: 

1 + tanh[sinh(u)] 
t= 2 ' 


( 9 ) 


To use the integral representation of Eq. (7) in practice, some algebraic operations need to be performed. Looking 
at the expression for the left passage probability in Eq. (1), one sees that for our purpose we only need to evaluate 
2 -F-| (a, 6 , c; z) with real parameters a = 1/2, b = k/ 4, and c = 3/2. The argument 2 = — cot (cf>) 2 is also real. To obtain 
parameter values that are suitable for the numerical integration of the integral representation, we apply the following 
combination of known results [37, 38] for 2 E-|: 

To evaluate 2 Ei(a, 6 , a; z), we use 


2 F 1 (a,b,a;z) = (1 - z) b . 

To evaluate 2 Ei(a, 6 , 6 ; z), we use 

2 F 1 (a,b,b-,z ) = (1 - z)~ a . 


( 10 ) 

( 11 ) 


Since we have a = 1/2 > 0, we apply the following transformation, based on the recursion relations for the Gauss 
hypergeometric function, to lower the value of a: 


CL — 1 — C 

2 Fi(a, 6 , c; z) = - -ry -— 2 Ei(a — 2, 6 , c; z) 


(a - l)(z- 1 ) 
c — 2 (a — 1 ) + (a — 1 — b)z 


2 F 1 (a- 1, b, c; z). 


( 12 ) 

(13) 


(a - l){z- 1 ) 

In case n < 4, we have b < 1, and apply the following recursion relation to increase the value of b in the integral: 


2 Ei (a,b,c;z) = ^ -— 2 Fi(a, b + 2 , c; z) 


b + 1 — c 
c — 2(6 + 1 ) + (6 + 1 — a)z 
b + 1 — c 

Using Eq. (15), we can end up with c < 6 ; in this case, we apply 


2 Ei(a, 6 + 1 , c; z). 


2 Fi(a, 6 , c; z) = —— C ^^ 2 -El (a, 6 , c + 2; z) 


c(c + 1)(1 - z) 

—c + ( 2 c + 1 — a — b)z 

c( 1 - z) 


2 Ei(a, 6 , c+ l;z), 


(14) 

(15) 

(16) 
(17) 


raising the value of c. 

The method for evaluating the left passage probability based on the integral representation has been implemented in 
C++ and is available online [42]. The class schramm_equation located in the header file schramm_equation.hpp takes 
the value of k in the constructor and provides a bracket operator which returns the numerical value of P K (c/). Applica¬ 
tion of the recursion relations and processing of special parameter values, as discussed above, is handled by recursive 
function calls in the header gauss_transf ormations .hpp. This, in turn, calls the class gauss Jiypergeometric which 
performs the numerical integration. 

Figure 3(a) shows the left-passage probability as function of />, for different values of kappa, as evaluated by the 
method discussed in this section. To test our numerical evaluation of the left passage probability, we compare the 
results obtained by our C++ program to results obtained using the Scipy (Python) [43] implementation of 2 F 1 , which 
uses the argument transformations described in Sec. II A. An example of this comparison is shown in Fig. 3(b) which 
plots the absolute difference between P 4 ((/>), as obtained from the presented method and from the Python library 
function, as function of <j>. 


III. CONCLUSION 

We discussed methods for the numerical evaluation of the Gauss hypergeometric function 2 Ei. While the presented 
options seem to be satisfactory for the parameter regime relevant for the evaluation of the Schramm-Loewner evolution 
(SLE) left passage probability, in general the numerical evaluation of 2 Ei is a complicated task and subject of ongoing 
research [32]. Hopefully, our work will facilitate numerical studies of SLE in the future. Given its importance for 
science, it would be useful if authoritative C++ libraries, such as the standard library [27], Boost [28], or GSL [35], 
provided implementations of 2 El that covered the relevant parameter regimes. 
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FIG. 3. (a) Left-passage probability P K (<(>) as function of 4>, evaluated using the method described in this work. For k, = 8, 

the curve is horizontal and it becomes more inclined with increasing k. The considered k values range from 1 to 8 as in Fig. 2. 
(b) Absolute value of the difference between P 4 as described in section IIB, and a control method using the Scipy (Python) 
2 F 1 implementation, as function of the angle tj). 
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